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Abstract 

A comparative simulation study of polymer brushes formed by grafting at a planar surface 
either flexible linear polymers (chain length Nl) or (non-catenated) ring polymers (chain length 
Nr = 2Nl) is presented. Two distinct off-lattice models are studied, one by Monte Carlo methods, 
the other by Molecular Dynamics, using a fast implementation on graphics processing units (GPUs). 
It is shown that the monomer density profiles p{z) in the z-direction perpendicular to the surface 
for rings and linear chains are practically identical, pr{2Nl,z) = Pl(Nl,z). The same applies to 
the pressure, exerted on a piston at hight z, as well. While the gyration radii components of rings 
and chains in z-direction coincide, too, and increase linearly with Nl, the transverse components 
differ, even with respect to their scaling properties: R g ^ y oc N^f 2 , Rgxy oc N® A . These properties 
are interpreted in terms of the statistical properties known for ring polymers in dense melts. 
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I. INTRODUCTION AND MOTIVATION 



Dense layers of macromolecules grafted by special chemical groups to non-adsorbing sub- 
strates have found enormous interest recently (for reviews see [EHO]) in view of numerous 
applications: colloid stabilization [7]; improvement of lubrication properties [2], also in a 
biological context [8] , creation of functional surfaces with switchable properties [9] , improve- 
ment of drug biocompatibility [ID], microffuidic devices for biomolecule separation [11], etc. 
The theoretical understanding of the (often unexpected) properties of polymer brushes, due 
to an interesting interplay of monomer- monomer and monomer-solvent interactions and the 
entropic forces resulting from the confinement of macromolecular conformations, has been 
a longstanding challenge as well (e.g. [TJ HJ |6j, fl2HT8] ) . 

To the best of our knowledge, all this research has focused on the grafting of linear 
polymers by a special end-group at one chain end; only as one exception the formation of 
"loop brushes" (where both chain ends are grafted to the substrate) has been considered [19] . 
In the latter case, although the chains are permanently entangled (forming a "catenated" 
network-like structure), the brush properties are found to differ at best marginally from 
those of brushes with free ends [19] . 

In recent years, there has been a great interest in the properties of polymer melts formed 
from non-catenated ring polymers [2DH4T)] . Both synthetic polymers (such as polyethylene 
or polystyrene [22]) and biopolymers such as DNA can be prepared as closed rings, and are 
of great interest as model system, to understand such diverse problems as polymer melt 
dynamics when the standard reptation mechanism (that needs chain ends [41] ) is eliminated 
[2D] EH EH], and the organization of DNA in the cell nucleus [371H2H35]. Note that most of 
the short genomes as well as plasmids are circular [16J H7J and also actin can self-assemble 
into rings [48] . Ring polymers under various kinds of confinement are also under discussion 
in this biological context (see, e.g., [4"D] |5"D]). 

It is a challenging problem to understand the conformation of ring polymers (both with 
respect to single collapsed rings [51] and confined rings jlDJ EDJ, E2] and rings in melts [SPI - 
ED]). One famous picture is that collapsed rings (as well as rings in a melt) are "crumpled 
globules" , where each subchain of the loop is condensed in itself, and the fractal dimension 

1/3 

of the object is df — 3. Hence the gyration radius scales like a compact globule, R g oc , 
Nr being the number of (effective) monomers in the ring polymer (hereafter, referred to 
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as ring length). However, when one hypothesizes that the structure of a ring polymer in 
a melt resembles a "lattice animal" (where every bond of the lattice is traversed by the 
polymer twice, in opposite directions) one predicts R g oc N^ 5 . The evidence from both 
experiment and simulations on this issue has been discussed controversially over decades, and 
only recently a resolution seems to emerge pEU E3> ESS HQ] : the ring polymers in melts behave 
like lattice animals for intermediate ring lengths but then crossover to crumpled globule- 
like structures for very long rings. Recall that both powers are substantially different from 
the law seen for melts of linear polymers, where a scaling like for gaussian chains occurs, 
R g oc N L |53j. We also stress that a (hypothetical) melt of ring polymers where the rings 
were free to cross each other, topological non-crossability constraints being "switched off" , 
would have trivial Gaussian statistics |53j as well. 

In the present paper, we ask what happens when ring polymers are grafted to a planar 
substrate by special groups at a single monomer. It is well known that for a polymer brush 
formed from linear chains at large enough grafting densities a g strongly stretched polymer 
configurations result while in the lateral direction the chains behave ideal [TJ [5] 

R$K*rN L ,R$ y <xa-^N]!* , (1) 

where Nl is the number of (effective) monomers in the linear polymer. The prefactor de- 
scribing the o" g -dependence can already be understood from the simple Alexander-de Gennes 
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[T2l [13] blob picture, describing a chain by a string of blobs of diameter d oc o g , each 
blob containing g = d 5 ^ 3 monomers, taking n = Nl/ g steps in the ^-direction perpendicular 
to the substrate while the typical excursion parallel to the substrate is of order y/nd. The 
self-avoiding walk statistics inside a blob, d oc g u with v = 3/5 is, of course, implicitly 
included in this argument. By comparing our simulation results for ring polymers with ring 
lengths Nr = 2 Nl of otherwise identical linear polymers forming a brush, at a g ^ = <?g,L/2 
(to keep the number of effective monomers strictly the same), we are able to clarify the effect 
of the topological interaction on the rings in this situation. We shall show indeed certain 
characteristic differences between ring brushes and linear brushes, which can be traced back 
to these topological effects. 
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FIG. 1: (a) Density profiles along the z-direction for rings and linear chains of different length (MD 
data). The density profiles scale linearly with the chain length for linear chains respectively half 
the length for rings. In the inset the curves collapse using this scaling relation, (b) Total monomer 
number density p(z), and end-monomers, p e (z) , respectively, middle- monomer density p m (z) from 
MC data. Here of = 0.0625 (rings) and Og = 0.125 (linear chains), (c) Pressure, exerted by the 
polymer brush of rings (circles), and linear chains (triangles) on a piston at distance z from the 
grafting plane - MC data for of = 0.0625 (rings) and Og = 0.125 (linear chains). Here Nr = 128, 
Nl = 64. The density profile is included for comparison. 

II. MODELS AND SIMULATION TECHNIQUES 

We always choose L x L grafting plane and periodic boundary conditions in x,y direc- 
tions. Chains are grafted at a grafting density <j g ^ = 0.125 for the linear chains with the 
range parameter a of the Weeks-Chandler-Andersen (WCA) potential that describes the 
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interaction between effective monomers in our MD simulations as unit of length 

U WCA (r)=Ae[(a/r) 12 -( ( x/r) 6 + ^, r < r c = 2 1 / 6 a . (2) 

r being the distance between any pair of monomers (bonded or nonbonded ones). Uwca{t > 
r c ) = 0, so the potential is purely repulsive. As is standard jUGS] bonded monomers interact 
with the finitely extensible nonlinear elastic (FENE) potential in addition, 

C/ F EN E (r) = -0.5^1n[l-(r/it; ) 2 ] , r < R , (3) 

where £/fene(?" > Ro) = 0. Parameters chosen are e = 1, temperature T = LOe/ks and 
k = 30, Ro = 1.5. The mass of the particles m = 1 as well, so the MD time unit is 
t = a I >Jm/e = 1 as well. MD runs were carried out on GPUs using the HooMD Blue 
r3574 code ji55J. Note that due to the use of GTX480 GPUs, a speedup factor of roughly 
70 in comparison to a run on a single i7 CPU core for the present application is obtained. 
We use a standard dissipative particle dynamics (DPD) thermostat [56] with parameters 
7 = 0.5, r cut = 1.25 • 2 1//6 , and integration time step St = 0.005. The grafting of the first 
monomer to the substrate site in the z = substrate plane is realized by the same FENE 
potential as in Eq. O). Note that the linear dimension L is adjusted such that we always 
have either M = 51200 or M = 102400 effective monomers in the system, for chain lengths 
N L = 16, 32, 64, 128, 256 and 512, respectively (recall N R = 2N L , a gfi = a g , L /2 throughout). 
For the longest chains, an equilibration time of r eq w 2.1 ■ 10 8 was needed. 

MC simulations were done using a somewhat different off-lattice model that we describe 
in the following. Of course, the motivation for choosing two different models is to confirm 
the universality of our results, which should not depend on irrelevant simulation details, but 
are rather generic. We employ a model of flexible polymer chains, each polymer consisting 
of Nl or Nr beads connected by anharmonic springs. These springs are defined by the 
FENE potential, Eq. @, albeit with r = I — £q. Thus, £ is the bond length, which can vary 
between £ min < £ < £ ma x, and has the equilibrium value £q = 0.7, while £ max = £q + Ro, 
£ m in = £o — Ro, and Ro = 0.3. With these parameters £ max = 1 is the unit of the length. A 
short-range Morse potential is used to describe non-bonded interactions, 

U M (r) = e M exp [-2(r - r min ] - 2exp [-x(r - r min )\, (4) 

with em = 1, r min = 0.8, a = 24. Since for this model the ^-temperature is fc^G ~ 0.62, in 
the present study we work at temperature /c^B = 1.0, and are thus in the "good solvent" 




FIG. 2: (a) Squared radius of gyration (MD) parallel to the grafting plane on double logarithmic 
scale Rg Xy = R 2 gx + R 2 gy against chain length Nl (linear chains) and Nr/2 (rings), respectivly. Inset 
- Scaling of the squared radius of gyration R 2 gz perpendicular to the grafting plane. The scaling 
prefactors differ by a factor of 4 between linear chains and rings, (b) The same from MC data for 
= 0.0625 (rings) and <jg = 0.125 (linear chains). 

regime. The Monte Carlo procedure consists of local displacements of the monomers only, 
which are accepted according to the usual Metropolis criterion. Typically, systems containing 
32768 monomers are simulated about 6 -¥■ 10 x 10 6 MCS after equilibration, whereby various 
quantities of interest are sampled. 
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III. RESULTS AND DISCUSSION 



One of the central static properties of polymer brushes is the density profile p(z) of 
monomers in the ^-direction perpendicular to the brush (Fig. fiji, b). A striking finding 
is that for both models, the density profiles of linear chains and rings are identical, within 
statistical errors (which are larger for the MC calculation, due to the much smaller statistical 
effort). The insert of Fig. [IJi shows the profiles rescaling z by Nl, to demonstrate that the 
data converge to the strong stretching limit. For small Nl there are oscillations for small 
z and a "finite size tail" for large z: these "corrections to scaling" become negligible in the 
limit of very long chains. However, the variation of p(z) near the brush height h (where 
p(z) vanishes for large Nl) is somewhat steeper than expected from the strong stretching 
limit of self-consistent field theory [HI [15], which would imply that p(z)/p(0) = 1 — (z/h) 2 : 
the latter result requires that the density inside of the blobs does not exceed the semi-dilute 
regime [53]. However, one should note that melt densities would be (in our units [51]) of 
order p « 1 for the MD model and p ~ 2 for the MC model. It is clear that our choice 
of parameters (leading to p(0) ~ 0.3, see Fig. [IJi) corresponds to the regime of a rather 
concentrated solution rather than a semidilute solution (remember that Eq. [2] accounts for 
the effects of a good solvent only implicitly rather than explicitly, of course). 

Fig. [T]d also includes the distribution of free chains ends (for the linear chains) and of 
the mid-monomers with index i = Nr/2 for the rings (which have Nr = 2Nl monomers, 
monomers being labeled consecutively along the ring, i — 1 being the grafted monomer). 
While according to the picture of the "Alexander brush" [12] Pe{z) oc 8{z — h), all end 
monomers need to be at the brush end, the actual distribution p e (z) shows that the ends (or 
mid-monomers in the case of rings, respectively) can be anywhere in the brush (although 
only few are near the grafting plane, since p e (z) oc z for z <C h [H]). It is interesting to 
notice that in p e (z) there are slight but systematic differences between linear chains and 
rings: for the latter the distributions are significantly more sharply peaked, and the tail 
towards large z is less pronounced. This is not surprising, of course, since the mid-monomer 
is bound by two strands rather than a single one. In Fig. [1J we show the pressure, exerted 
by a brush or ring (circles), or linear polymers (triangles) on a piston at height z above the 
grafting plane. Evidently, the pressure closely follows the monomer density profile at the 
brush tail, whereby one harldy detects any difference between rings and linear chains. 
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A very interesting behavior results when we examine the linear dimensions of individual 
chains, however (Fig. [2^i, [2]d). While the z-components of the mean square gyration radii in 
z-directions (R 2 gz ) for linear chains and rings coincide precisely, the transverse components 
differ significantly, indicating a different power law, 

{R] xy )<xN (linear chains), (R^ xy ) oc N - 8 (rings) . (5) 

While the first of these equations is expected (Eq. (jl|), the second is not. Surprisingly, 
our data coincide with the theoretical prediction for lattice animals iy = 2/5) [26] over a 
full decade in N, for both models. For very small N, deviations occur, as expected, but 
indications for a crossover to the "crumpled globule" exponent [51] iy = 1/3) are not seen, 
unlike the case of bulk polymer melts [3T)1 - PIU] . Of course, while it is conceivable that another 
crossover could occur for much larger ring lengths, it is not really clear that one should expect 
this: after all, the components of the grafted chains in the brush are strongly stretched, 
(R gz ) oc N 2 , while in a ring polymer melt all components R gx , R 2 gy and R gz scale in the same 
way. Due to this strong anisotropy, it remains unclear if the transverse components R 2 gxy 
of the chains in a ring polymer brush should exhibit the same scaling as chains in a ring 
polymer melt. The fact that the mean-square gyration radii in ^-direction, (R 2 Z ), for linear 
chains and rings (under equivalent conditions, N R = 2N L , = &g/2) coincide, while the 
transverse components are so different in both cases, should not be mistaken as a proof that 
^-components and transverse components are decoupled, however. When one compares the 
full probability distributions P(R 2 Z ) for rings and linear chains, one does see characteristic 
differences: for very small R gz , the distribution in the case of rings has much more weight 
than it has for linear chains (Fig. [3]). Of course, this finding does not contradict Figs. [2^, |2}o, 
since the region of very small R 2 makes only negligibly small contributions to (R 2 qz )- 



IV. CONCLUDING SUMMARY 



In this letter, we consider a new class of polymer brushes, formed by densely grafted ring 
polymers to a flat substrate surface. We demonstrate that the transverse configuration of 
such grafted rings differs substantially from corresponding linear chains (Fig. [2j although 
collective properties (such as the average density profiles of the monomers (Figs. [lp,,b) or 
the pressure exerted by the brush on a compressing wall ((it)) basically do not differ at all. 
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FIG. 3: Normalized probability distributions of Rg Z for rings and linear chains of different length 
N. The shoulders in the probability distribution of the open chains towards lower values of Rg Z 
indicate that this distribution is non-gaussian - there is a higher probability of finding compact 
configurations of linear chains as compared to rings of double length. Inset: Normalized probability 
distribution of the squared radius of gyration in x-direction parallel to the grafting plane R 2 gx . 



An intriguing question is the interpretation of our findings in terms of geometric properties 
of chain versus ring conformations (Fig [4] shows some illustrative snapshots). In the corre- 
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sponding problem of non-catenated rings forming dense melts, one associates a R g ~ N R 
variation with lattice animal like conformations (while for very large ring-lengths Nr one 
rather finds (39] crumpled globules, leading to R g ~ N R ). Obviously, the projections of 
monomer positions into the xy-plane for rings (Fig. [1]) are fairly compact, while the corre- 
sponding projections for linear chains are random walk-like. 

An intriguing question is what phenomena are affected by these different conformations 
of rings (rather than linear chains) forming a polymer brush. Presumably, for instance, 
the dynamics of free (linear) chains from a solution or melt flowing past a ring polymer 
brush will be very different. The interaction with nanoparticle inclusions, that involves the 
deformation of chains (or rings) surrounding the nanoparticle would be rather different. All 
such questions remain to be explored also. 

Of course, it would be very nice if guidance from analytical theory were available to 
understand the differences between ring polymer brushes and ordinary brushes formed from 
linear polymers, as we have found in our studies. However, none of the well-established 
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FIG. 4: Configurations from Molecular dynamics simulations to illustrate the conformational dif- 
ferences between loops and open chains: 3 out of 400 loops of ring length Nr = 256 (left) and 3 out 
of 800 open chains of chain length Nl = 128 (right) are displayed. The rings and the open chains 
are grafted on a 80 x 80 grafting plane. The overall number of monomers is 102400 per system. 

approaches (self-consistent field theory [T4T1T6] , density functional-type theories [31 157] , etc) 
can take the topological differences between catenated and non-catenated rings into account: 
none of the theories could hence predict any of the differences that we have found here. Also 
experimental studies of this problem (ideally one migth consider neutron scattering from 
deuterated rings in a protonated ring polymer brush) are a great challenge. 
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